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ABSTRACT 

We consider the evolution of a warped disc around a Kerr black hole, under conditions 
such that the warp propagates in a wavelike manner. This occurs when the dimension- 
less effective viscosity, a, that damps the warp is less than the characteristic angular 
semi-thickness, H/R, of the disc. We adopt linearized equations that are valid for 
warps of sufficiently small amplitude in a Newtonian disc, but also account for the 
apsidal and nodal precession that occur in the Kerr metric. Through analytical and 
time-dependent studies, we confirm the results of Demianski & Ivanov, and of Ivanov 
& Illarionov, that such a disc takes on a characteristic warped shape. The inner part of 
the disc is not necessarily aligned with the equator of the hole, even in the presence of 
dissipation. We draw attention to the fact that this might have important implications 
for the directionality of jets emanating from discs around rotating black holes. 
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1 INTRODUCTION 

Accretion discs are found around black holes in the centres of 
active galaxies and also in galactic X-ray binary stars. These 
systems may form in such a way that the angular momen- 
tum vectors of the disc and of the black hole are not parallel. 
As a particle in a tilted orbit around a rotating black hole 
undergoes Lense-Thirring precession at a rate dependent on 
the radius of the orbit, so a tilted disc experiences a gravit- 
omagnetic torque that tends to twist and warp the disc. If 
the disc is able to communicate a warping disturbance in a 
diffusive manner as a result of its pressure and viscosity, it 
may adopt a characteristic warped shape in which the plane 
of the disc undergoes a smooth transition from one plane to 
another in the vicinity of a certain radius (Bardeen & Pet- 
terson 1975; Kumar & Pringle 1985; Scheuer & Feiler 1996). 
At large radius the disc is essentially flat and its plane is de- 
termined by the total angular momentum vector of the disc. 
At small radius the disc is essentially flat and lies in the 
equatorial plane of the black hole. Over a longer time-scale 
the disc and hole tend towards mutual alignment (Scheuer 
& Feiler 1996; Natarajan & Pringle 1998). 

However, there are circumstances in which a warping 
disturbance propagates in a wavelike, rather than diffusive, 
manner. This occurs in a Keplerian disc when the dimen- 
sionless viscosity parameter a (Shakura & Sunyaev 1973) 
is smaller than the angular semi-thickness H/R of the disc 
(Papaloizou & Lin 1995; Wijers & Pringle 1999; see the dis- 
cussion by Pringle 1999); in a non- Keplerian disc the prop- 
agation is also wavelike, although it is dispersive (Ogilvie 
1999). The evolution of a disc around a rotating black hole 



under these conditions has received less attention. On the 
basis of analytical calculations, Demianski & Ivanov (1997) 
and Ivanov & Illarionov (1997) found that the disc can adopt 
a steady warped shape in which the tilt angle is an oscilla- 
tory function of radius. More recently, Nelson & Papaloizou 
(2000) conducted smoothed particle hydrodynamic simula- 
tions of accretion discs around rotating black holes under 
a variety of assumptions. They did not find the tilt oscilla- 
tions, but found that the inner part of the disc was aligned 
with the equatorial plane of the hole. 

In this paper we re-examine the shape of a warped disc 
around a Kerr black hole, under conditions such that the 
warp propagates in a wavelike manner. We explore the rea- 
sons for the existence of a steady wavelike solution, and dis- 
cuss the expected wavelength, amplitude and phase of the 
oscillations in disc tilt. Using time-dependent calculations 
of the linearized equations, we study the process by which 
such a steady solution is established. Finally, we discuss the 
potentially important implications of the shape of the disc 
for the directionality of jets emanating from discs around 
black holes. 



2 WAVELIKE PROPAGATION OF WARPS 

The propagation of warping disturbances of small amplitude 
in a nearly Keplerian and nearly inviscid disc has been dis- 
cussed by Papaloizou & Lin (1995), Masset & Tagger (1996), 
Demianski & Ivanov (1997), Pringle (1999) and Lubow & 
Ogilvie (2000). These authors differ in respect of the nota- 
tions and assumptions they adopt, but all agree that bend- 



2 S. H. Lubow et al. 



ing waves propagate approximately non-dispersively under 
these conditions, with a wave speed that is somewhat less 
than the average sound speed of the disc. 

We consider disturbances of a steady disc of surface den- 
sity S(-R) and vertically integrated pressure P(R) defined by 



E = / pdz 



P = / pdz, 



(1) 



where (R,<j),z) are cylindrical polar coordinates. Let O(R), 
k(R) and O z (R) be the orbital angular velocity, the epicyclic 
frequency and the vertical oscillation frequency associated 
with circular orbits at radius R from the central object. After 
an integration by parts, and using the vertical hydrostatic 
equilibrium of the disc, we find 



(-3) 



dz — Q, 



pz 2 dz = 2 z T,H 2 , 



where H(R) is an effective density scale-height. In the case of 
a vertically isothermal disc, H is equal to the usual Gaussian 
scale- height. 

A simple way to describe a warped disc is as a continu- 
ous set of tilted circular rings, and to define a unit tilt vector 
l(R,t) normal to the plane of the ring of radius R at time 
t. When the disc is only slightly tilted out of the a;y-plane, 
l z ~ 1. As shown by Lubow & Ogilvie (2000), the equa- 
tions governing warping disturbances of small amplitude in 
a nearly Keplerian and nearly inviscid disc may be written 
in the form 



dG 

dt 



dt R OR 
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20 



lxG + aOG 



PR 3 dl 
4 dR' 



(3) 



(4) 



Equation ([]) is derived from the vertical component of the 
equation of motion, and can be understood as expressing the 
conservation of horizontal angular momentum: 2irG(R,t) is 
the horizontal internal torque in the disc, and T(R, t) is the 
horizontal external torque, per unit area, acting on the disc. 
Equation (Q) is derived from the horizontal components of 
the equation of motion, and determines the internal torque 
G. The internal torque is mediated by shearing epicyclic 
motions, driven by the radial pressure gradients that are set 
up when a stratified disc is warped. These motions undergo 
apsidal precession when k ^ O, and also viscous decay, here 
described by a dimensionless viscosity parameter a. 

The parameter a appears in equation (Q) only to quan- 
tify the rate, aO, at which the shearing epicyclic motions as- 
sociated with the warp decay. Although this term is derived 
from a viscous model, numerical simulations of the inter- 
action of shearing epicyclic motions with three-dimensional 
magnetorotational turbulence (Torkelsson et al. 2000), and 
an analytical closure model of such turbulence (Ogilvie 
2002), both suggest that small-amplitude shearing epicyclic 
motions decay exponentially in time. The measured rate is 
in fact through to be somewhat smaller than aO if a is cali- 
brated in the usual way with reference to the i?<j!>-component 
of the turbulent stress. 

It is often more convenient to adopt a complex repre- 

Then 



sentation in which W = l x + il y and G = G x + iG y 
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Figure 1. Dispersion relation for bending waves in the case of 
prograde nodal and apsidal precession with Q^/fi = 0.95 and 
k/Q = 0.9, according to equation (R) (solid lines) and equation 
m (dashed lines). 
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Note that the external torque is related to the non-sphericity 
of the potential, and therefore to the difference between O 2 
and O z . In addition, W(R,t) represents the tilt of the disc at 
each radius in the sense that W(R,t) = /3(R,t) exp[i'y(R,t)], 
where /3 is the amplitude of the local tilt and 7 is the azimuth 
(cf. Pringle 1996). 

These equations are known to be valid when the quan- 
tities W and G are sufficiently small, and vary on a length- 
scale long compared to H and on a time-scale long compared 
to O^ 1 . In deriving them it is also assumed that the quan- 
tities |1 - 2 z /0 2 \, |1 - k 2 /0 2 \ and a are of order H/R or 
smaller. 

In the absence of viscosity, the dispersion relation asso- 
ciated with this system, for wavelike solutions with a rapidly 
varying phase factor exp J (iu dt — ik dR), is 



O 2 - Oj 
20 



O 



20 



= -k 2 H 2 2 z . 
4 



(7) 



This shows that, in an exactly Keplerian disc, the warp prop- 
agates as a non-dispersive wave with wave speed HO z /2, 
which is equal to half the isothermal sound speed in the 
case of a vertically isothermal disc. The presence of nodal 
and/or apsidal precession introduces dispersion, and forbids 
wave propagation for frequencies between (fi 2 — fi^)/2^ and 
(O 2 -k 2 )/20 (Fig. 1). 

It is interesting to compare this dispersion relation with 
that derived by Lubow & Pringle (1993), which describes 
the full set of wave modes of a vertically isothermal disc. 
Ogilvie & Lubow (1999) slightly extended that work to allow 
for the possibility that O z 7^ O, and found that, for a disc 
undergoing isothermal perturbations (7 = 1), 
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k 2 H 2 



(8) 



where & = oj — mft, m is the azimuthal mode number and 
n is the vertical mode number. For the case m = 1, n = 0, 
which corresponds to a tilt or warp, we find 



ft 2 - ft 2 + u/ 



2fi 

ifc 2 .ff 2 ft 2 
4 



n 2 - K 2 + 0/ 

2n 



(9) 



Unlike equation (Q), this dispersion relation is valid when 
the quantities |1 — Q 2 /Q 2 \, |1 — k 2 /Q 2 \ and kH are of order 
unity. The two equations clearly agree well in the limit of 
low frequency, |w/f2| <C 1, and agree exactly when ui — 
(see Fig. 1). 

The case u — is of some interest, as it relates to 
warped discs that are independent of time. When both nodal 
and apsidal precession are present, and are either both pro- 
grade or both retrograde, the dispersion relation for ui = 
indicates that the spatial structure of the warp is oscillatory 
in character, its radial wavenumber being given by 



k 2 H 2 



(q 2 - n 2 z )(n 2 - k 2 ) 

ft 2 ft 2 ' 



(10) 



An oscillatory solution of this kind may be interpreted as a 
bending wave having zero phase velocity but non-zero group 
velocity. 

Conversely, when the nodal and apsidal precession are 
in opposite senses, a steady warp is spatially evanescent. 
This situation is qualitatively similar to the case when the 
viscosity is large and the warp satisfies a diffusion equation. 
It applies when the orbital, epicyclic and vertical frequencies 
are derived from an axisymmetric Newtonian gravitational 
potential that satisfies Laplace's equation in the mid-plane 
of the disc, 



2JT + D.i 



0. 



(11) 



since this implies 

sgn(ft 2 - ft 2 ) = -sgn(ft 2 - k 2 ). (12) 
However, this need not hold in the metric of a black hole. 



3 THE STEADY SHAPE OF A WARPED DISC 
AROUND A BLACK HOLE 

The orbital, epicyclic and vertical frequencies around a Kerr 
black hole satisfy (Kato 1990) 



n 2 



( GM 



) (r 3/2 + a) 



(13) 



ft 2 

n 2 - rt'i 



= 6r 



8ar z ^ 2 + 3a 2 r 



-3/2 



o 2 -2 

6a r , 



ft 2 ( 15 ) 

where — 1 < a < 1 is the dimensionless angular momentum 
parameter of the hole (a < implying a retrograde disc) 
and r is the radius in units of GM/c 2 . For a prograde disc 
(a > 0) we have the ordering k < Q z < ft at any radius, 
corresponding to a situation as illustrated in Fig. 1. 
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Figure 2. Ratio of the radial wavelength of the steady warp 
to the effective density scale-height, plotted against the distance 
from the black hole, expressed in units of GM/c 2 . The four curves 
relate to Kerr black holes with a = 0.1 (top), a = 0.2, a = 0.4 
and a = 1 (bottom). The inner radius plotted corresponds to the 
marginally stable circular orbit in each case. 



Fig. 2 shows the radial wavelength, A = 2n/k, associ- 
ated with the steady wavelike shape of a prograde warped 
disc around a Kerr black hole. To a first approximation, i.e. 
considering only the first terms in equations wo 
have 



H 



r 5/4 



(6o)Va 



(16) 



The wavelength is shorter for more rapidly rotating black 
holes, and decreases as the marginally stable circular orbit 
is approached, but is always at least a few times H. This 
suggests that these wavelike solutions are physically realiz- 
able and may be described with some confidence by theories 
of long-wavelength bending waves. 

A more accurate description of the steady wavelike 
shape can be obtained by setting the time-derivatives to 
zero in equations (^) and (^|). This results in a second-order 
ordinary differential equation for the tilt variable W(R), 



_d_ 

dR 



pr 3 q 



Q 2 - /t 2 + 2iaft 2 



AW 
dR 



+ T,R 3 (D 2 -n 2 z )W = 0.(17) 



If we take the first-order approximations from equations 
©-©> we have 



T.H 2 



3i — 1 + ia 



AW 
AR 



+ 8ar~ 3/2 EW = 0. 



(18) 



Under the further assumptions that E tx RT 1 ^ 2 and H tx 
R, this equation is equivalent to that studied by Ivanov & 
Illarionov (1997), who first identified the steady wavelike 
solutions. 

A WKB treatment of equation (^) in the absence of 
viscosity also indicates that the radial variation of the am- 
plitude of the oscillations in W is proportional to 



4 S. H. Lubow et al. 



tl 2 



n 2 



-1/4 



(19) 



(20) 



and therefore, to a first approximation, proportional to 
(Eff) 1 / 2 ' 

The amplitude may either increase or decrease as one ap- 
proaches the inner radius, depending on the properties of 
the disc, but is likely to do so slowly. 

More specifically, if E oc RT 2a and HJR = er h_1 with 
e = constant and h > a > — | , equation dlq) in the absence 
of viscosity has a solution in Bessel functions, 



W = x v [C 1 JAx) + C 2 Y v {x)], 



where 

h — a 



h + 



/24a\ 1 / 2 r-^+W 
{ — ) h+ i ■ 



(21) 



(22) 



At large R (x <C 1) the second Bessel function is dominant 
and the tilt tends to a constant value, 



W^W OB = -—V(v)C: 



(23) 



At small R (x 3> 1) the Bessel functions are oscillatory and 

2 \ 1 / 2 



W 



\nx / 
+C2 sin 1 v 



Ci cos I x — 



~2 



(-f-l) 



(24) 



Note that the amplitude of W oc x"~ 1/2 oc R 1/S (XH)- 1/2 
as quoted above. In the special case considered by Ivanov 



& Illarionov (1997), a 
amplitude of W oc R~ 
when v = |, since then 



and h 



1, so v 



and the 



1/8 



A different special case occurs 



W = 



2 x 1/2 



(I) 



(Ci sin x — C2 cos x) 



(25) 



exactly. 

While the constant C2 is fixed by the tilt at large radius, 
C\ is determined by the inner boundary condition. If this 
is such that the torque vanishes at the inner radius R = 
-Rin (usually at the marginally stable circular orbit), then 
dW/dR — there. In the special case v = this determines 
the solution as 

cos(:Ti n — x) 



(26) 



COS X ln 

provided that cos X[ n 7^ 0. This solution can be understood 
as a standing wave consisting of a superposition of ingoing 
and outgoing bending waves (having zero phase velocity but 
non-zero group velocity). The condition of vanishing torque 
at the inner radius reflects the ingoing bending waves. Note 
that in this case the warp consists of pure tilt, without twist- 
ing (i.e. 7 = constant). 

The steady shape of the warp can be interpreted in 
terms of a resonant wave- launching process (e.g. Goldreich & 
Tremaine 1979). However, in the present case, the resonance 
is at R = 00 because the forcing, caused by the steady disc 
tilt at large radius, occurs at zero frequency. In addition, the 
region of maximum forcing, close to the black hole, is far re- 
moved from the resonance (cf. Lubow & Ogilvie 2000, where 
the resonance is at R = and the forcing is concentrated 



in the outer part of the disc). As a consequence, the forcing 
effectively occurs at an intermediate radius, typically where 
the wave executes its outermost wavelength (kR ~ x ~ 1). 
The wave is launched as a trailing (k > 0) inwardly propa- 
gating wave, having group velocity less than or of order the 
sound speed (see also the discussion by Nelson & Papaloizou 
2000). 



4 TIME-DEPENDENT CALCULATIONS 

We now demonstrate how these equations can be applied 
to the time-dependent behaviour of a disc in the neigh- 
bourhood of a rotating black hole. The simple application 
we have in mind here is the Newtonian approximation to 
Lense-Thirring precession around a rotating black hole. In 
this case, we can write approximately (Demianski & Ivanov 
1997; Nelson and Papaloizou 2000; Section 3 above) 



V = 
and 



o 2 



2fi 2 



3^ 
2 R 



n 2 z -n 2 
2n 2 



a ( R B 



3/2 



(27) 



(28) 



where _R S = 2GM/c? is the Schwarzschild radius. 

In the numerical calculations we choose to solve for the 
quantities A(R,t) and D(R,t), which define the radial ve- 
locity perturbation and the enthalpy perturbation according 
to 



u R 



Aze 14 



Dze 4 



(29) 



where the real part is understood. These are simply related 
to W and G by W = -D*/Rtt 2 z and G = ^EH 2 R 2 nA*. To 
be explicit, we use the equations (|H]) and ^) to solve for A 
and D. Thus, 



OA 
~dt 
and 

dD 
~dt 



D 

£i 
2 



1 dD 

2 7m 



+ irjQA - aVLA, 



d -(T,H 2 R 1/2 A) 



(30) 



(31) 



E// 2 R 1 /2 QR v 

We consider a disc which is inclined to the spin axis of 
the hole at large radius. For ease of comparison with the an- 
alytic results we consider the particular case in which steady 
shape of the disc is described by Bessel functions of order 5 , 
and so can be written in terms of trigonometric functions. 
Thus we assume the disc with inner radius at R = R ln 
to have E = T, in (R/ R in )~ 1/2 , H = H in (R/ R in ) 3/4 , Q = 
^in(R/ Rin)~ 3 ^ 2 , and thus (in the vertically isothermal case) 
a sound speed c s = c Si i n (_R/_Ri n ) _3//4 . If we let 

7 7 = ?7in ( J R/i? in )' 1 , C = (in(R/Ry s/2 , 
then the steady solution of the equations, with a = 0, is 



(32) 



D^P^R 
and 

i/3oc.A 



A 



277inf2in-Rii 



COs[A/i?in] 



■(R/Rin) 



cos[A/_Ri n 



(33) 



(34) 



where 
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,2 _ 4r; in C in _R in 
= HI 



(35) 



and Poo is the tilt angle of the disc as R — > oo. This solution 
is chosen to ensure that there is zero torque (^4 = 0) at the 
inner boundary. 



4.1 Numerical results 

The numerical method is to take A to be defined at the N 
logarithmically distributed grid points (typically N = fOOf 
or 2002), and to take D to be defined at the half-grid points. 
We take A and D to be defined both at full time points and 
at the intermediate half-time points. A is updated from: (a) 
D at the intermediate time-point using straightforward nu- 
merical differencing for the radial derivative, (b) A at the 
intermediate time-point for the precession term (leapfrog), 
and (c) A at the previous full time-point for the dissipative 
term (Euler method). Similarly D is updated from: (a) A at 
the intermediate time-point using straightforward numeri- 
cal differencing for the radial derivative, and (b) D at the 
intermediate time-point for the precession term (leapfrog). 
Boundary conditions are required only for A, and we take 
A — at the boundaries, which corresponds to a zero torque 
condition at the inner and outer disc radii. 

For numerical convenience, we take units in which Ri n — 
1, Poo = 1, and Qi n = f . We then choose the other param- 
eters to produce a Newtonian approximation to a rotating 
black hole for which the angular momentum parameter is 
a — | (4 — \/T0) w 0.5585 and for which the marginally 
stable circular orbit is at R ln — 2R B . This implies that 

Vin 



-§ and On 



1(4- VW) 



-0.1396. We then 
choose H in /Ri n = [(4 - v / T0)/87r 2 ] 1/2 ^ 0.1030 to ensure 
that A = 27ri?i n . This gives a steady solution (with viscosity 
zero, i.e. a — 0) of the form: 



D = RT 2 cos[27r(l - R- 1 )]. 



(36) 



This solution has the property that the tilt angle 
f3(R) = \DR 2 \ is unity at R = 1, at R = 2, and as R -> oo, 
and is zero at R — 4 and at R — 4/3. As a check, this so- 
lution was tested numerically by using it as initial input on 
a grid from 7?i n = 1 to R ou t = 30, using N = 1001 logarith- 
mically spaced grid points. No evolution was detected. 

We then investigated time-dependent behaviour, by 
considering the effect of taking a disc which is initially 
aligned with the spin axis of the hole, and tilting the outer 
parts. Thus as an initial condition we take DR 2 = for 
R < 18, DR 2 = i(l + sin[7r(i? - 20)/4]) for 18 < R < 22, 
and DR 2 = 1 for R > 22. To consider a reasonable physi- 
cal situation we add a small amount of dissipation and take 
a = 0.05. We take A — initially, which satisfies the inner 
and outer boundary conditions and also corresponds to the 
correct solution for a disc with uniform tilt around a Kep- 
lerian point mass. Thus the initial condition corresponds to 
a discontinuity in the disc tilt at radius R = 20 smoothed 
over a distance of AR — 4. We integrated the equations for 
a time of t = 20, 000 using a grid with N = 2002 extending 
from R = 1 to R = 900. 

In Fig. 3 we show the initial evolution up to time 
t = 4, 000. As expected the discontinuity in the tilt initially 
propagates inwards and outwards in a wavelike fashion. The 
outward propagation continues throughout the calculation. 
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Figure 3. The tilt of the disc is shown as a function of radius at 
11 equally spaced times t = 0, 400, 800, . . . , 4, 000. The tilt at the 
final time, t = 4, 000, shown in bold, shows that the steady state 
solution has already been established at radii R < 8. 
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Figure 4. The tilt of the disc is shown as a function of radius 
is shown for the inner disc regions at 11 equally spaced times 
t = 0, 400, 800, . . . , 4, 000. The tilt at the final time, t = 4, 000, 
shown in bold, shows that the steady state solution has already 
been established at radii R < 8. 
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Figure 5. The tilt of the disc is shown for the inner disc regions 
at 11 equally spaced times t = 0, 2000, 4000, . . . , 20, 000. The tilt 
at the final time, t = 20, 000, shown in bold, shows that the steady 
state solution has already been established at radii R < 65. 



The inward propagation of the tilt interacts with the pre- 
cessional effects at the hole and by the time t = 4, 000 the 
steady state shape of the disc close to the hole is essentially 
established (Fig. 4). It takes the form of a steady warp, to- 
gether with a slight twist (caused by the small amount of 
viscosity) as predicted by Ivanov & Illarionov (1997). The 
interaction of the inwardly propagating warp with the in- 
ner disc regions generates a reflected wave which propa- 
gates outwards. Thus towards the end of the computation 
(t = 20,000), shown in Fig. 5, at which time the initial 
outwardly propagating warp wave is approaching the outer 
edge of the grid, followed closely by the reflection of the ini- 
tially inwardly propagating warp wave, the steady warped 
disc solution has been established over about half the grid. 



5 DISCUSSION 

We have considered the time-dependent evolution of a 
warped centrifugally supported disc in the low-viscosity 
limit (a<H/R) in which warp propagation is wavelike 
rather than diffusive. We note that this is relevant to discs 
in which the viscosity is low, such as protoplanetary discs, 
and/or in which the disc thickness is large, such as non- 
radiative or high-luminosity discs around black holes. We 
make use of the linearized perturbation equations, which re- 
quire that the tilt angle of the warp be sufficiently small. 

As an illustration we apply the resulting equations to 
a Newtonian approximation of a moderately thin (H/R ~ 
0.1), low viscosity (a — 0.05<-ff/f?) disc around a Kerr 
black hole with angular momentum parameter a « 0.5. The 
disc in this case is subject to both apsidal and nodal pre- 
cession, and the evolution equations we derive are similar 
to those derived by Dcmianski & Ivanov (1997) in the limit 



of low viscosity. Ivanov & Illarionov (1997) demonstrated 
that these equations give rise to a steady state in which the 
disc takes up a characteristic warped shape with the pre- 
cession induced by the potential being balanced by wave- 
like stresses in the disc. We confirm this result, and we use 
our time-dependent equations to demonstrate how such a 
steady state is set up, starting from the initial condition of 
a disc whose inner parts are aligned with the spin axis of 
the hole, and whose outer parts are tilted at a fixed incli- 
nation. We show that, once the outer tilt has propagated 
(in a wavelike fashion) to the central regions of the disc, the 
steady-state configuration becomes rapidly established. The 
adjustment of disc shape at the inner radii is propagated 
outwards through the disc at the usual local wave propaga- 
tion speed. 

Another possible approach to describing the steady 
shape of a warped disc around a black hole is based on the 
work of Ogilvie (1999), who developed a non-linear theory 
of bending waves in accretion discs. The derivation is based 
on certain ordering assumptions that eliminate some of the 
time-derivatives from the problem. In the time-dependent 
case, this theory is applicable to the generic, 'non-resonant' 
case and not to the special case of discs that are both 
nearly Keplerian and nearly inviscid. However, in the case 
of steady warps, the theory of Ogilvie (1999) is applicable to 
all regimes and has the advantages of being fully non-linear, 
and of including the effect of the accretion flow. We do not 
pursue this idea here, but note that for an inviscid disc with 
a steady warp, the theory implies that 

ldG +T, (37) 
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(38) 



where 7 is the adiabatic exponent and \ip\ = R\dl/dR\ is 
a dimensionless measure of the amplitude of the warp. In 
the linear regime these equations agree with equations (^|) 
and (^) with the time-derivatives set to zero. The theory 
indicates that non-linear effects become noticeable when 
\RdW/dR\ 2 is comparable to (SI 2 - k 2 )/D. 2 . This can be 
expected to occur in the inner part of the disc unless the tilt 
at large radius is sufficiently small, i.e. /3oo <C 1. The effect 
of the non-linearity is likely to be to increase the wavelength 
of the warp and restrict its amplitude. 

Other authors have cautioned against taking seriously 
the possibility of the disc warp having a substantial ampli- 
tude in the centre of the disc. Nelson & Papaloizou (2000) 
did not find such warps in their non-linear SPH compu- 
tations and attribute this failure to a presumption that 
"non-linear effects lead to the damping of short wavelength 
[warps], and thus cause the alignment of inner disc regions 
in which the tilt amplitude would otherwise change rapidly 
on small length-scales". Although this may be the case, we 
note that for the case we consider, with H/R ~ 0.1, the ra- 
dial length-scales of the warps are generally at least of order 
the local radius, R, (Fig. 2) and are not small compared to 
the local disc scale H . We note further that once the radial 
length-scale of the warp variation becomes short, the disper- 
sive nature of the warp waves is likely to become significant 
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and may delay the onset of non-linearity of the waves. In 
addition, in their original paper, Ivanov & Illarionov (1997) 
point out that although the formal computations predict a 
disc warp at inner radii, such predictions might be nullified 
if the warp itself gives rise to turbulence within the disc, 
and so a high local effective viscosity. They suggest that the 
oscillatory vertical shear within the disc caused by the warp 
(Papaloizou & Pringle 1983) might become shear unstable. 
Gammie, Goodman & Ogilvie (2000) have investigated this 
possibility, and show that instability occurs through a para- 
metric effect and is likely to set in when \A\ >30aQ. 

We have not discussed the interesting question of the 
net torque between the disc and the black hole, or the time- 
scale for mutual alignment, under conditions such that the 
warp propagates in a wavelike manner. The torque may be 
considered to be exerted in launching the steady train of in- 
wardly propagating bending waves, which carries a certain 
flux of angular momentum. Unlike the case of resonantly 
launched waves in problems of discs subject to periodic tidal 
forcing (Goldreich & Tremaine 1979), the torque cannot be 
simply expressed in terms of the properties of the disc in the 
neighbourhood of a certain radius. Furthermore, if the vis- 
cosity is small enough that the waves reach the inner radius 
and reflect from it, the torque that would cause a mutual 
alignment is partially or completely cancelled. For example, 
in the case v = i leading to the solution ( pt^ ) in the ab- 
sence of viscosity, the integrated horizontal torque is given 
in term of an integral over the disc of the relevant compo- 
nents of the local torque T(R, t). Writing T = T x +iT y , the 
total horizontal torque is 



2tt / TRdR 

'Rin 



COS x ln 



IWoo tan Xix, 



C0s(xin — x) dx 



(39) 



Rather than causing a mutual alignment, this torque causes 
a slow mutual precession of the disc and black hole. The 
direction and magnitude of the precession depend sensitively 
on the parameters of the disc. 

Finally, we note that if the outer disc and the black 
hole are misaligned by more than 90° so that the disc can 
be considered retrograde, the steady wavelike solution is re- 
placed by an evanescent solution. This occurs because the 
nodal and apsidal precession are in opposite senses (see 
equations 0, and |l|). In that case the steady shape of 
the disc may be expected to be qualitatively similar to that 
envisaged by Bardeen & Petterson (1975), with the inner 
disc aligned with the equator of the black hole but rotating 
in the opposite sense. We illustrate this in Fig. 6. 



q 0.6 




Figure 6. The result of a calculation identical to that of Fig. 3, 
but in which the sign of the nodal precession f is reversed. This 
corresponds to the physical situation in which the disc and hole 
rotate in opposite senses. 



without changing the tilt of the outer disc. These results 
contrast with the usual finding (e.g. Nelson & Papaloizou 
2000) and/or assumption (e.g. Natarajan & Pringle 1998) 
that the inner regions of the disc align with the equator 
of the hole. In the discussion above, we have noted above 
that confirmation of these results awaits a proper calcula- 
tion using full general relativity, as well as an assessment of 
possible non-linear, dispersive and parametric effects. Nev- 
ertheless, it is evident that the results presented here could 
have important implications for the directions in which jets 
might emanate from accreting spinning black holes. Indeed, 
if the region responsible for direction of jet collimation is 
at several radii from the hole, which is likely to be the case 
for relativistic jets, then the Newtonian approximations ap- 
plied above may be adequate to confirm the effect, even if 
the very inner regions of the disc are indeed aligned by the 
fully relativistic effects close to the hole. A lack of corre- 
lation between the inner and outer disc tilts could provide 
one explanation of the finding by Kinney et al. (2000) that 
the directions of jets from low luminosity AGN appear to be 
uncorrelated with the disc plane of the host spiral galaxies. 



6 CONCLUSION 

The result that, in low- viscosity discs around a Kerr black 
hole, the inner parts of the disc are not necessarily aligned 
with the black hole, as found by Ivanov & Illarionov (1997), 
is a general one. Indeed, depending on the disc properties, 
it is possible for the inner disc to be tilted at a greater an- 
gle to the hole than the outer parts of the disc. In addition, 
because the inner disc shape depends sensitively on the ra- 
dial dependence of disc properties (such as surface density 
and disc thickness), a change (for example) in the accretion 
rate can give rise to a change in the inner disc warp, even 
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